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It is argued that some phase-transitions observed in models of non-equihbrium wetting 
phenomena are related to contact processes with long-range interactions. This is investigated 
by introducing a model where the activation rate of a site at the edge of an inactive island of 
length £ is 1 + ai~°' . Mean-field analysis and numerical simulations indicate that for cr > 1 
the transition is continuous and belongs to the universality class of directed percolation, 
while for < cr < 1, the transition becomes first order. This criterion is then applied to 
discuss critical properties of various models of non-equilibrium wetting. 

I. INTRODUCTION 

Many recent theoretical studies have shown that the growth process of a solid phase on a sub- 
strate can undergo a variety of non-equilibrium transitions. They are analogous to equilibrium 
wetting phenomena in which liquid boundary-layers exhibit critical behavior in the vicinity of the 
liquid-gas coexistence line. Such growth processes can be effectively modelled by defining a suitable 
evolution rule for the profile (interface), corresponding to the boundary of the solid layer. Several 
growth models which do not obey detailed balance and evolve towards stationary non-equilibrium 
states have been studied in the past. In many cases, by varying a control parameter, they exhibit 
a transition from a regime where the solid phase remains pinned to the substrate to a regime 



where an unbounded growth sets in 



This depinning of an interface may 



be considered as non-equilibrium wetting, in analogy with its equilibrium counterpart. In some 
models the character of the depinning transition changes, depending on the dynamical rates which 
control the interface evolution. Numerical studies in 1 -|- 1 dimensions allowed to identify both first- 
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FIG. 1: Schematic configuration of an interface bound to a substrate. Within the DP framework bound 
sites may be considered as active while islands of depinned sites are inactive. 

and second-order phase-transitions, those of the latter type falling into two universality classes: 
directed percolation (DP) and multiplicative noise (MN) Since the transition always oc- 
curs towards an irreversibly depinned phase, it is quite natural to draw an analogy with dynamical 
processes characterized by an absorbing state. 

Recently, it has also been shown that such a depinning transition can be related to the synchro- 
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12l |. In this latter framework, 



nization phase transition in spatially extended chaotic systems 
two different replicae of the same dynamical system are coupled one to each other, either deter- 
ministically 12] or by the addition of the same realization of a spatiotemporal stochastic noise. 
Upon increasing the coupling parameter, the system undergoes a non-equilibrium phase transition 
between an unsynchronized phase and a completely synchronized one, characterized by a vanish- 
ing local difference between the two (initially) different replica. Since two completely synchronized 
replica of the same system are bound to remain identical one with respect to each other at all future 
times, we can conclude that the completely synchronized phase is absorbing for system dynamics, 
thus playng the role of the irreversibly depinned phase in non-equilibrium wetting processes. Inter- 



estingly, both numerical analyses of coupled map lattices systems |lll . Il2l | and analytical arguments 



12 



13l | predict the synchronization transition to exhibit a critical behavior, which belongs either 
to the MN or to the DP universality class. 

A schematic typical configuration of an interface bound to a substrate is depicted in Fig. ^ It 
is composed of detached domains separated by pinned segments. The dynamics of the interface is 
such that each detached domain may either shrink or expand from its edges; new domains may be 
created by the unbinding process of bound sites, and two or more detached domains may merge 
into a single larger one. In principle, a segment internal to a detached domain may bind back to 
the substrate. However, in some physical conditions such processes are virtually suppressed. This 
occurs when the unbound interface moves on the average away from the substrate, while it is held 
bound to the substrate by some short range attractive interaction. In fact, the farther is a segment 
from the edge of a domain, the larger is the height of the inferface and, accordingly, the more 
unlikely the possibility to bind back to the substrate. In this case, the dynamics of the interface 
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may very well be described by a contact process in which the active sites correspond to those 
bound to the substrate. The resulting depinning transition is thus expected to belong to the DP 
universality class. The observation of DP depinning transition in some models llJI confirms 
the validity of the above arguments. In some cases a possible crossover to a first order transition 
has been suggested, when the attractive interaction between the interface and the substrate is 
increased. While in certain regions of the phase-diagram the existence of a first order transition 
has been clearly demonstrated, in other regions it has only been tentatively suggested, based on 

, . Due to the strict connection between wetting and 
complete synchronization non-equilibrium phase transitions, a more accurate understanding of this 
part of the wetting phase diagram is highly desirable not only from the theoretical point of view, 
but also in view of possible experimental realizations of the MN and/or the DP universality classes. 
In our opinion, in fact, the synchronization transition in spatially extended systems represents the 
most promising framework in which one can look for MN critical properties, if not for DP ones, 
which seem to require a highly non linear local dynamics governing a spatially extended system 
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1^ . It is worth stressing that both these classes are still eluding a clear experimental evidence, 
possibly due to the presence of quenched disorder in experimental realizations and, in the case of 
DP phase transitions, to seemingly unavoidable small fiuctuations of the (supposedly) absorbing 
state Interesting candidates for experiments include the photosensitive Belousov-Zhabotinsky 
(BZ) reaction, which is known to exhibit complex spatiotemporal dynamics and semiconductor 
lasers with time-delayed optical feedback. In Ref. it has been shown that such delayed 

systems can be interpreted in terms of a suitable spatiotemporal dynamics, where the effective 
system size is given by the ratio between the delay time and the typical fast timescale of the 
system. Notice that, in principle, the use of semiconductor lasers with time-delayed feedback 
allows for obtaining rapidly a large data set, which is basically free of quenched noise effects. 
These features provide the possibility of a proper statistical description of the synchronization 
transition. Finally, the role played by small fiuctuations near the completely synchronized phase 
and their exponential suppression in the case of synchronization DP has been discussed in Ref. |l^ . 



In the present paper we introduce a framework within which the crossover from DP to first 
order transition may be examined. This framework is then applied to two previously introduced 
models of non-equilibrium wetting. 

Preliminary simulations performed in the pinned phase close to a seemingly first-order wetting 
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transition have revealed that the activation rate at the border of a depinned island depends on 
the island length. This suggests that the dynamics of a fluctuating interface leads to an effective 
interaction between the sites at the island-boundaries. If this interaction is long range, it can 
affect the dynamics of large islands and, in principle, provide a mechanism for a first order wetting 
transition. In this paper we consider a generalization of the contact process where the activation 
rate of sites at the boundary of an inactive domain decays algebraically with the domain length. 
In particular, the activation rate for an island of length i is assumed to take the form A(l + a/i"), 
where A, a and a are positive constants. We find that, depending on the power o", the model 
exhibits either a continuous DP-like (for o" > 1) or a first-order (for cr < 1) transition. We 
then examine the possible emergence of such effective long range interactions in wetting models. 
Numerical determination of the power a of the effective interaction thus allows inferring the order 
of the phase transition also in these models. 

The paper is organized as follows: In Section 2 we introduce a generalized contact process 
in 1 + 1 dimensions that includes an activation rate which decays algebraically with the size of 
inactive domains. The mean-field solution of the model predicts a first-order transition, when 
the interaction decays slowly enough with the inactive domain size, while faster decay rates yield 
a DP behavior. In Section 3 we present the results of detailed numerical simulations, which are 
in a very good agreement with the mean field predictions. The bridge between this generalized 
contact process and wetting models is discussed in Section 4. There we consider two models, which 
have been introduced previously for studying non-equilibrium wetting. The first is a solid-on-solid 
model P, y, 0], whose phase-diagram contains a line of wetting transitions. While the first order 
nature of the transition has clearly been established in a part of the phase-diagram, the nature of 
the transition in another region proved to be more difficult to analyze. Our numerical study seems 
to suggest that inside the latter region, upon changing a suitable control parameter, the effective 
activation rate at the boundary of inactive domains can exhibit both a fast (cr > 1) and a slow 
(fj < 1) power law decay. However further analysis based on scaling arguments shows that the 
slow decay of the effective activation rate (and thus the first order behavior) is only a finite-size 
effect, albeit particularly robust. This result suggests that in this entire region the transition is 
asymptotically continuous and of DP nature. We then analyze a second wetting model (single- 
step- with- wall) 0, where previous studies have indicated a continuous DP-like wetting transition 
for strong attractive interaction between the substrate and the interface. We provide numerical 
evidence that the effective activation rate governing the dynamics of inactive islands corresponds 
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to the case a > 1. This is indeed consistent with the DP nature of the transition. The main 
results are summarized in Section 5, which contains also some remarks and comments on future 
perspectives. 



II. A CONTACT PROCESS WITH LONG RANGE INTERACTIONS 



In this section we introduce a lattice model of a contact process in 1 + 1 dimensions with long 
range interactions which shows a crossover between a continuous DP and a first-order phase- 
transition. We consider a periodic lattice of length L, where the state variable Si at site i is either 
"active", Si = 1, or "inactive", Si = 0. The dynamics evolves by random-sequential updates, i.e. 
at each time step a lattice site i is chosen at random. If the selected site is either active or next 
neighbor of an active site, it is updated according to the following rules 



1^0 with rate 1 (1) 
oh [10^] ^ O^-^l [110^-^] with rate A(£) = A(l + a/r) , 

where 0^ is a shorthand notation for an inactive island of size £. Finally, as for usual contact 
processes, inactive sites that are not adjacent to active ones cannot be activated, thus guaranteeing 
that the inactive state is absorbing. The constants A and a are both non-negative: the case 
a = corresponds to the usual short range contact process, which exhibits a DP transition at 
Ac = 1.64892(8) For a > the power law decay with £ of the shrinking rate of inactive islands 
introduces effective long-range interactions. Within the mean field approximation one finds that 
the transition is continuous for a > 1 but becomes first-order for < cr < 1. To demonstrate this 
point let p{t) be the average density of active sites at time t. In the thermodynamic limit, L oo, 
the mean field dynamics of p reads 

I = -p + Xp^ + ^)(l - pY = (A - l)p - Xp^ + Xap^ ■ (2) 

e=i 1=1 

For cj > 1 the sum on the r.h.s. of Eq. @ is finite in the limit p — > and its contribution amounts to 

a renormalization of the coefficient of the p^ term. Accordingly, the mean field equation describes 

a standard DP process with short range interaction, thus recovering, for sufficiently small values of 

a, the continuous nature of the phase-transition. For large a, the coefficient of the p^ term becomes 

positive and the mean-field approximation predicts a first-order transition. However, studies of 
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similar models indicate that the mean-field prediction of first-order phase-transitions associated 
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with the change of sign of the term is unreliable 

For < cj < 1, the leading contribution arising from the sum in the rightmost side of Eq. @ 
can be captured by replacing it with an integral over 

dtp={X-l)p-Xp^ + Xap^ ^-/^de, (3) 

which, to leading order in p, reduces to 

dtp = (A - l)p + aAr(l - a)p^+'' - Xp^ . (4) 

Here T{x) is the standard Gamma function. The leading nonlinear term in this equation involves 
a non-integer power, as a consequence of the long range nature of the interactions. Since its 
coefficient is positive, Eq. ^ cannot admit fixed point solutions for arbitrarily small densities, so 
that the transition to the absorbing state is first order. This result is expected to hold even for 
(7 = 1, where the leading singular term in the equation is — In p. 



III. NUMERICAL RESULTS 



The mean field calculation, that is expected to hold above the upper critical dimension of 
directed percolation dc = 4, indicates that for a > 1 the transition is second-order with the 
critical exponents of DP, while for o" < 1 it turns into a first-order transition. The crossover 
from a continuous to a discontinuous transition is therefore predicted to take place at ctc = 1. In 
what follows we investigate the validity of this prediction by simulating the generalized contact 
process defined in Eq. In order to obtain independent checks, we performed two different 
kinds of numerical analyses, i.e. measuring the scaling properties of suitable observables starting 
from i) a fully active state and ii) a single active site (epidemic spreading). Although we do 
not expect the overall scenario to depend on the parameter a (see Eq. ), the accuracy of the 
numerical simulations actually does depend. For small a, discontinuities in the order parameter are 
correspondingly small, while for large a, it is necessary to consider very large lattice sizes to reach 
the asymptotic-scaling regime. All the numerical results reported in this paper have been obtained 
for a = 2, which represents a good compromise. We find (Tc = 1.0 it 0.1, which suggests that the 
mean-field analysis is quite accurate even in 1 + 1 dimensions. However, we cannot exclude the 
possibility that the precise threshold value in such a low dimensional case slightly deviates from 
1. This has already been observed in directed percolation with long-range infections through Levy 
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Ai 



0.5 



0.. 



0.9 



1.1 



1.2 



1.5 



1.2850(5) 1.3650(5) 1.4091(1) 1.4708(3) 

1.2870(4) 1.3395(3) 1.3635(5) 1.4093(3) 1.4280(5) 1.4711(3) 



TABLE I: Estimates of the critical points for a = 2. They have been obtained starting from a fuUy active 
initial state (A{) and by spreading analysis of localized initial conditions (A^). Finite size corrections are 
responsible of the small differences. 



flights, where a smaU deviation from the mean-field prediction was found for the critical value of 
the control parameter (i.e., the exponent of the Levy distribution) Ifll. l2(]|. 



A. Analysis of the stationary active state 

We first discuss the evolution of Monte Carlo simulations that start from a fully active state. 
In order to distinguish between first-order and second-order transitions we determine the density 
of active sites and the size-distribution of inactive islands. Let us first summarize the expected 
results for each of the two quantities. Since the continuous transition should be DP, the average 
density of active sites, p{t), measured at criticality, should decay as 

Pit)-t~-' • (5) 

On the other hand, off-criticality, in the active phase of an infinite system, p{t) saturates to a 
stationary value {p)t (where {■)t denotes time average), which scales with the distance from the 
critical point Ac as 

(p)i~|A-A/ , (6) 

where /? = 0.276486(8) 2l|- 

These two critical exponents are connected by the so called temporal 
exponent zyy, so that = jd/uw = 0.159464(6) I2l']. 

Moreover, in DP the size-distribution P{tj of inactive islands decays algebraically as 



P(£) ~ ^-(^-/^/^'i) (7) 

with a cutoff at the spatial correlation length ~ |A — \c\~^^, where the exponent 2 — (3/v± = 
-1.747. . . follows from {i) ~ l/p ~ i'^/'"^ . In fact, the average length of inactive islands 

poo 

{£) = / P{e)U£ (8) 

^0 
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FIG. 2: Power law critical behavior of the generalized contact process Q for cr = 1.1 and a — 1.5, compared 
to the expected DP behavior (dashed lines): a) decay in time of the density of active sites p(t); b) stationary 
density of active sites as a function of the distance from criticality. Circles correspond to the case a — 1.5, 
while squares to a = 1.1. Both graphs are plotted in a doubly logarithmic scale. 

diverges at criticality due to scale invariance, so that P{i) ~ with 1 < a < 2 for i ^ Since 
the average size of inactive islands is proportional to the inverse of the density of active sites, one 
has 1/p ~ {1} ~ ^1^'"^, which implies a = 2- f3/u± = 1.747 . . . Q- 

Conversely, for a first-order phase-transition p{t) is not expected to exhibit any critical behavior 
associated with a diverging correlation length at the transition point. Instead, the saturated order 
parameter {p)t exhibits a discontinuity at the transition. Since in this case the active phase cannot 
display any coarsening properties, the average length (jSJ of inactive islands should be finite, i.e. 
P{i) should decay faster than As no hysteretic behavior can be observed in non-equilibrium 

processes with an absorbing state, the distribution P{£) turns out to be the most effective indicator 
of a first-order transition. 

In order to reduce as much as possible finite-size effects, we considered very large systems of 
size L = 2^^ with periodic boundary conditions (additionally we have also averaged over a few 
different realizations). The best estimates of the critical point Ac, are reported in the first row of 
TableHl For a = 1.1 and a = 1.5 we identified a critical scaling region (see Fig. I^J, where both the 
exponents 5 and P are in agreement with the best numerical estimates for DP in 1+1 dimensions. 
Moreover, we computed P{i) in the active phase close to the critical point by sampling spatial 
configurations at periodic time intervals and counting all the inactive regions of size i. As shown 
in Fig. |SJ -P(^) is characterized by a power law decay slower than 1/^^ (prior to the unavoidable 
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FIG. 3: Doubly logarithmic graph of the (im-normalized) size distributions of inactive islands P{£) as 
a function of size £. From top to bottom the solid lines corresponds to ct = 1.5, a = 1.1, a = 0.9 and 
(7 = 0.5. The dot-dashed line marks the power-law decay expected in the case of a DP phase-transition. 
The long dashed line decays as discriminating in 1-1-1 dimensions between first order and continuous 

phase-transitions. 

exponential cutoff due to finite-size effects), consistently with the prediction for DP (see Eq. 0). 

On the other hand, below a = 1 a scaling region at the transition point could not be identified. 
The saturated density of active sites {p)t shows a finite discontinuity and P{i) decays faster than 
(see Fig. ^ . Accordingly, this analysis provides evidence that the transition is discontinuous. 



B. Spreading from a single seed 

A further verification of the mean-field analysis has been obtained by simulating model (0, 
starting from a single active site at the origin 25]. In this type of simulations, the relevant 
variables are the survival probability Ps{t), the number N{t) of active sites (averaged over all runs) 
and the mean square spreading R'^{t) of the active region. At the critical point of a phase-transition 
towards an absorbing state, these quantities are known to scale as 

p,(^)~^-^ N{t)r^t\ R^{t)^t'^i'. (9) 

In the special case of a DP phase-transition, the exponents 5, ry, z can be expressed in terms of the 
standard exponents /?, u^, v\\ as 



(10) 
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FIG. 4: Typical cluster grown from a single seed for a = 0.5 at criticality, simulated up to 5000 Monte 
Carlo sweeps. Notice that large inactive islands are suppressed and the growing cluster can be regarded as 
effectively compact. 





b T\ z 


DP 


0.159464(6) 0.313686(8) 1.580745(10) 


Glauber 


1/2 2 



TABLE II: Exponents in seed simulations for DP and zero-temperature Glauber dynamics in 1+1 dimen- 
sions. 

Their actual values are reported in the first hue of Table ^ A scaling behaviour of the type 
described by Eg. is expected to arise also in systems exhibiting a first-order transition in 
1-|-1 dimensions although the relations between spreading and stationary exponents assume 
a more general form than Eq. (|injl In these cases, the critical dynamics follows from the 
marginal relative stability of two coexisting phases: (i) the absorbing state itself and (ii) ^ phase 
characterized by a suitable finite density of active sites. For instance, an active site evolves into 
a droplet of phase (ii) embedded into a sea of phase (i). At criticality, the boundary, having no 
preferential velocity, diffuses like an unbiased random walk. A clear such example of such a type is 
the one-dimensional Glauber-Ising model at zero temperature, where the density of active sites in 
phase (%%) is maximal, i.e. equal to 1. In models such as the natural existence of small inactive 
islands (see Fig. |1} can seemingly make the separation into two phases questionable. However, for 
sufficiently small o"- values, the long-range interactions are strong enough to suppress the formation 
of large inactive islands. In these circumstances, one would expect the long-term dynamics to 
be controlled by the evolution of the boundaries. In particular, the scaling exponents of contact 
processes exhibiting a first-order phase transition in 1-1-1 dimensions are expected to be the same 
as in the Glauber-Ising model at zero temperature, whose values are reported in the second line 
of Table ^ Notice that r] = means that the density remains finite at criticality. 

Numerical simulations of the spreading dynamics confirm our expectations. Independently of 
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FIG. 5: Time evolution of the average number of active sites close to the phase transition. In the left panel 
the average number of active sites N{t) is plotted as a function of time for a — 2, a — 0.5,0.8,0.9 and for 
different values of A, close to Ac- All curves indicate that N{t) converges to a constant after a transient 
time which increases with a. In the right panel the time-rescaled number of active sites N{t)t^''' (with 
f] = 0.313686) is showed for a — 2, a = 1.1,1.2,1.5, and for different values of A, close to Ac- In order 
to avoid a messy overlap, the curves have been shifted vertically by an arbitrary value. Their asymptotic 
behavior indicates good agreement with the expected DP critical scaling. The estimates of Ac are reported 
in Table m Both graphs are plotted in a doubly logarithmic scale. 

the results discussed in the previous section, we have first estimated Ac by measuring the average 
number of active sites N{t) for different values of A and then looking for the value of the control 
parameter that minimized the curvature of N{t) at long times (see Fig. By performing simu- 
lations up to tmax = 3 X 10^ Monte Carlo sweeps and by averaging over 10^ realizations, we have 
obtained fairly accurate estimates, which are listed in the second row of Tabled The differences 
with the values reported in the first row provide an indirect estimate of the magnitude of finite-size 
corrections. 

We have also investigated the asymptotic behavior of Ps{t) and N(t) at criticality (see Fig. IB]). 
For large times, both quantities exhibit a power-law behavior: for a > 1 (a < 1) the growth rates 
are consistent with a DP transition (Glauber- Ising dynamics). As expected, the closer is a to 1, 
the longer is the time needed to reach the scaling regime. 

In conclusion, our simulations indicate that the crossover from a continuous to a discontinuous 
transition takes place in our model between a = 0.9 and a = 1.1. This result is consistent with 
the predictions of the mean field arguments discussed in Section 2. For the sake of completeness. 
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FIG. 6: The survival probability Ps{t) (left panel) and the average number of active sites N{t) (right 
panel) at criticality for cr = 0.5,0.8,0.9, 1.1, 1.2, 1.5 (showed from bottom to top) in a doubly logarithmic 
representation. The predicted asymptotic slopes of directed percolation (DP) and Glauber-Ising (GI) are 
indicated as bold lines. 

we mention the results of simulations made in the marginal case a = 1. In such a case, there 
are indications of a still first-order transition, although the Glauber-Ising exponents are not yet 
recovered on the accessible time scales. 



IV. NON-EQUILIBRIUM WETTING AS A CONTACT PROCESS WITH LONG 

RANGE INTERACTIONS 



In this section we investigate to what extent the behavior of non-equilibrium wetting processes 
can be interpreted as a contact process with long-range interactions of the form and whether the 
previous results can be used as a criterion for distinguishing first-order from DP-like continuous 
transitions. To this end we consider two previously introduced wetting models, studying them 
within the above derived framework. In practice, we numerically estimate the effective activation 
rates at the boundary of detached islands and show that they are indeed of the form The 
results obtained in the previous section can help to discern the nature of the wetting transition in 
the two models. 
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A. Restricted solid on solid wetting model 

The first system we consider is a restricted solid-on-solid (RSOS) model Q, IJ]. It is defined 
on a one dimensional lattice with periodic boundary conditions: at each lattice site i the height 
variable hi can take any nonnegative integer value such that — = 0, 1. A hard attractive 
substrate is located at zero height preventing hi from becoming negative. The interface evolves by 
random sequential updates controlled by three real parameters q, qo and p. At each move, a site 
i is randomly selected and, provided the above constraints are fulfilled, one of the following three 
processes is carried out (see Ref. 

- Particles are deposited with rate qo at the bottom layer and with rate q at higher layers. 

- Particles evaporate from the edges of a terrace with rate 1. 

- Particles evaporate from the middle of a plateau with rate p. 

The phase-diagram of the model is shown in Fig. [3 If qo is large enough (e.g. qo ^ q), the model 
exhibits a line of continuous phase-transitions, which belongs to the MN universality class. Its 
critical behavior can be described by a Kardar-Parisi-Zhang (KPZ) equation in a potential 
V{h) representing the interaction between the interface and the substrate 

dV 

h = DV^h + u{Vhf - — + (. (11) 

Here h = h{x,t) indicates the height of the interface on the substrate, while C = Ci^it) is the 
noise term (5-correlated in space and time, {C{x,t) ^{x' ,t')) ~ 5{x — x')5{t — t'). The coefficient 
1/ of the nonlinear term is positive for p > 1 and negative for p < 1. Detailed balance holds only 
for p = 1. This special case can be solved exactly (28|: its critical properties are described by an 
Edwards- Wilkinson (EW) equation [2^ {u = 0) equipped with the interface-substrate interaction 
potential. 

Decreasing qo amounts to increasing the attractive interaction between interface and substrate. 
It has been observed that, below a certain threshold, the continuous transition may turn into a 
discontinuous one ,3]. For instance, at p = 1 it was shown analytically that for < 2/3 the 
transition becomes first-order, while the transition point remains located at q = 1. Numerical 
simulations provide very clear evidence that this scenario extends to the non-equilibrium case 
p > 1: the phase-transition line is still independent of qo (see Fig. ITJ, while for go smaller than 
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FIG. 7: Phase diagram of the RSOS model for qo > Qoip) (left panel) and go < Qoip) (right panel). The 
phase-transition takes place along the full line, while the dashed line in the right panel marks the lower 
border of the phase-coexistence region (see text). 

a threshold value Qq^p), the transition becomes first-order. In this case, the interface dynamics 
is quite different from the one described by our model. In fact, direct inspection of the interface 
dynamics in the active phase close to criticality shows that there is a non-negligible probability for 
the interface to return in contact with the substrate not only at the domain boundaries of inactive 
domains, but also inside these domains. In DP jargon, this amounts to saying that an inactive site 
may become active even without being in contact with an active site. This excludes any direct 
relation with the model of a generalized contact process with long-range interactions introduced 
in Section 2. 

Conversely, for < p < 1, the dynamics appears to be strongly related to that of our model. For 
sufficiently small qq (i.e. qo smaller then a threshold Qq^p)), a region in the {p, q) plane arises, where 
the pinned and the unbound phases coexist (see Fig. [7jp). In this region, an unbound interface 
moves away from the substrate and never binds back. On the other hand, a bound interface remains 
bound for macroscopically long times and, in the thermodynamic limit, will never detach from the 
substrate. Only if the growth rate q is increased beyond a new critical value, a depinning transition 
takes place at the upper border line of the phase-coexistence region (see Fig.jTJa). It is the nature 
of this transition that we want to investigate here. 

As explained in P, the stability of the pinned (active) phase in the coexistence region is 
ensured by the negative sign of the nonlinear coefficient v: when a large detached island forms, 
it grows quickly, acquiring a triangular shape with a given slope and eventually shrinks from the 
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FIG. 8: RSOS model - decay of the activation rate (see text) as a function of size of depinned islands. 
Numerical simulation have been performed close to the special point p = 0. Circles refer to Tq = (p = 
0, qo = 0.2, q — 0.755), while squares to Ti ~ {p ~ 0.01, qo = 0.2, q ~ 0.76). The graph is plotted in a doubly 
logarithmic scale. For the sake of clarity the data are shifted vertically by an arbitrary value. 

edges with constant velocity. Because of the triangular shape of the interface, the probabihty of 
returning to the substrate at some point far from the edges of the island is exponentiahy suppressed. 
Therefore, it is reasonable to conjecture that for p < 1 the model belongs to the class of contact 
processes. 

Numerical simulations of the RSOS model dynamics (not shown here) seem to indicate that 
also in this case the depinning transition may become discontinuous for qq < q^ip)- For p = and, 
correspondingly, go < 9o('-*) ~ 0.399, we know that it is of DP type 3]. On the other hand, the DP 
scaling regime becomes transparent only after a transient time that, for go = 0.35 is on the order 
of 10"^ units. Although the crossover becomes practically unobservable for yet smaller go-values, 
there are compelling reasons to believe that this regime is eventually attained P]. 

In order to shed some light about the order of the transition for < p < 1, we have tested 
whether effective long-range effects spontaneously emerge as a result of the RSOS microscopic 
dynamics. In practice, we have measured the effective activation rate A(£) = Na{i)/Niy{£) at the 
border of depinned islands of size i {Niy(l) is the number of times a depinned site at the border of an 
inactive island of size i has been selected in the stationary regime and Na{i) is the number of times 
the selected site is immediately pinned). Data has been obtained by averaging over time, space 
and different realizations (typically, 100) for large lattices (L = 10^) and close to the transition 
line. 
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FIG. 9: RSOS model - decay of the activation rate (see text) as a function of size of depinned islands. 
Numerical simulation have been performed at the critical points T2 = (p = 0, qo ~ 0.4, q 0.70) (squares) 
and Ta = (p = 0.2, qq = 0.3, q = 0.745) (circles). The graph is plotted in a doubly logarithmic scale. Data at 
p — 0.3 show a power law decay with an exponent a = 0.69(2), while qo = 0.4 data decays with an exponent 
cr = 0.89(2). 

We have first analyzed the nature of the transition close to p = 0. In Fig. E] we display the 
results obtained at the transition points Tq = {p = 0, qo = 0.2, q = 0.755) and Ti = {p = 0.01, qo = 
0.2, q = 0.76). The activation rate A(^) is reported after subtracting its estimated asymptotic 
value Aoo = lim£__^oo ^(^)- In both cases, it is found that it converges towards Aqo faster than 
although no precise estimate of the scaling rate can be obtained. Altogether, these results suggest 
that DP critical properties persist also for small values of p and qo < qoip). On the other hand, 
for p = 0.2 and qo < qo{0.2) = 0.515... our findings are suggestive of a first-order phase transition 
(and thus, in agreement with previous findings |3|)- In fact, from Fig.lHl we see that in both points 
T2 = {p = 0.2, qo = 0.4, g = 0.70), and T3 = (p = 0.2, go = 0.3, g = 0.745) the activation rate 
is found to converge slower than l/l, although the actual value of the exponent a < 1 appears 
to depend on the parameter. According to the criterion introduced in Sec. 2, these results are 
therefore compatible with a discontinuous phase-transition. 

As a result, we can attribute the seemingly first-order nature of the transition to the existence of 
effective long-range interactions. Nevertheless, it remains to be proved whether the slow dependence 
of the activation rate on the window size is just a finite-size effect or holds fot arbitrarily large 
distances. Here below we present an argument supporting the former hypothesis. 

In order to clarify this point we move progressively away from the equilibrium case. For p = 1, 
it is known that the transition is discontinuous [J] (for small enough qo) and that the dynamics of 
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the free interface is asymptotically described by the EW equation. When p is lowered below 1, the 
only relevant difference that is expected to occur is a crossover in the free interface dynamics from 
and EW to a KPZ regime above some critical length ic 

It is thus natural to conjecture 

that as long as the dynamics of the bound interface remains insensitive to such differences, the 
scenario should not change. Here below we argue that this occurs for system sizes smaller than 
some length Lc that can be extremely large. 

The best way to characterize the above mentioned crossover is by monitoring the width w{t) = 
{{{h — {h))'^))^^'^ of an initially flat (and free) KPZ interface. In fact, after an initial growth as t^^'^ 
(in agreement with EW equation), at some time tc, w{t) crosses over towards a behavior of the 
type t^^^ and eventually saturates because of the finite length of the system. In other words, for 
times smaller than tc the interface behaves in the same way as in an equilibrium regime. Since the 
stationary profile of an interface is a diffusive random walk, one can safely assume that 



where L is the interface length. By extending the above relation to finite times, one can interpret it 
as the definition of the effective scale L{t) that is resolved at time t. For instance, ic = {w{tc)/k)'^ is 
the minimal length of a free interface that allows observing a crossover to the KPZ scaling behavior. 

Within the context of a bounded interface, this implies that deviations from equilibrium are 
observable only in those depinned islands of length i > ic- Accordingly, the problem of determining 
the minimal length to observe deviations from the equilibrium scenario amounts to estimating 
the probability for a suitably large depinned islands to arise. At equilibrium, the theoretical 
analysis developed in Q has revealed that when the transition is discontinuous (the scenario we 
are interested in), the interface is exponentially localized at the substrate, i.e. the probability to 
find large values of the height h scales as P{h) ~ exp{—h/ho). It is quite plausible to assume, and 
we have indeed numerically verified, that as long as deviations from equilibrium are not detectable, 
the exponential decay survives also for p < 1. Since we have also seen that depinned islands have an 
approximately triangular shape, this means that also island lengths are exponentially distributed 
P(i) ~ exp{—i/io), where £o is proportional to ho, the proportionality constant being related to 
the slope of such islands. 

As a result, the probability that at least a given island reaches the size ic is proportional to 
exp(—£c/io)- In a large but finite system of size L, this may happen independently at different 
places. Hence the first large island would appear in a typical time r ~ ic/L exp(^c/^o)- Accordingly, 




(12) 
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the minimum system size guaranteeing that such islands are observed with nonnegUgible probabihty 
and dominate the wetting dynamics is 

L, = 4 exp(4/4) • (13) 

Therefore, as long as the (effective) size remains smaller than Lc, a seemingly first-order transition 
is observed. Beyond Lc, in the fully nonlinear regime, several theoretical and numerical studies of 
different models j^,^^ suggest that a DP behavior sets in. As a result, we expect in particular that 
the distribution of depinned islands crosses over from an exponential to a power-law distribution. 
However, given the exponential dependence of Lc on ^c, it may happen that the crossover length 
is so large that the asymptotic regime is practically unobservable. Note also that diverges 
exponentially as p — > 1, thus approaching the equilibrium point. 

Numerical simulations performed at the transition point Ti = {p = 0.01, g = 0.76, go = 0.2) 
indicate that k = 0.179(1) and w{tc) = 2.9(2), yielding the estimate ~ 260. On the other 
hand, direct computations of the island sizes indicate that £o ~ 0(10^), yielding Lc ~ 3.5 x 10^ 
and altogether confirming that the crossover towards DP can be observed, as we actually do. 
Moreover, at the transition point T2 = {p = 0.2, g = 0.70, go = 0.4), we find k = 0.204(1) and 
w{tc) = 5.7(2), while £0 is almost unchanged, yielding the estimates lc ~ 780 and Lc ~ 2 x 10^. 
Consistently, no indication of the crossover has been observed up to times on the order O(IO^) and 
system sizes of length 0(10^). 



B. Single Step with Wall model 

The second model we have tested is the so-called "single-step- with- wa 
a variant of the well known Single-Step model introduced in Refs. 



30|, 



1" (SSW) model 



3ll |. Here, the growing 



interface is described by a set of integer heights hi at site i of a one-dimensional lattice of length 
L with periodic boundary conditions. In this model, the "continuity" restriction — = 1 
plays the role the RSOS constraint. An upward-moving wall is located at some integer height 
hw{t), below the interface. It moves with velocity v^, thus pushing the interface which cannot be 
overtaken by the wall. Moreover, in analogy with the RSOS model, the interface is also attracted 
by the wall. The model evolves by random-sequential dynamics, i.e., at each time step dt = 1/L, 
a site i is chosen at random. If the interface has a local minimum at site i (namely hi < /ij±i), 
the height hi is increased by two units with probability 1 if /ij > /i^, or with probability (1 — g) 
whenever the interface is pinned to the wall {hi = h^)- Since < g < 1, an effective attractive 
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FIG. 10: Updating rule of the SSW model. The full line represents the interface, while the shaded area 
represents the wall. Dashed segments indicate interface growth occurring in randomly chosen local minima 
(see A and B in upper panel) and in all sites located below the wall after it has been shifted upwards by 
one unit (see C in lower 
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FIG. 11: SSW model - Double logarithmic plot of the activation rate A(^) as a function of the size £ of 
depinned islands. Numerical simulations have been performed at criticality: circles corresponds to q = 
0.7, Vw = 0.4297 while squares to q ~ 0.8, v^u = 0.3489. Dashed lines mark our best fits with Eq. 
rendering respectively the estimates a — 1.20(5) and a = 1.30(5). 

force is introduced between the wall and the interface. After = L/vw time steps, the wall is 
moved upward by one unit, while the height of all interfacial sites that would be overtaken by the 
wall is increased by two units (e.g., see Fig. Il()|) . On the basis of these microscopic update rules, 
one can easily infer that the velocity of the free interface is 1/2 in the thermodynamic limit. 

The phase-diagram of the SSW model is controlled by two parameters: v^ — l /2, i.e. the relative 
velocity of the wall to the free- interface, and g, i.e. the "stickiness" of the wall. By decreasing Vyj, a 
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phase-transition from a pinned to a depinned phase is observed. Numerical analysis and analytical 



arguments show that when q < q* = 0.444. . . such a transition takes place at = 1/2, i.e. 
when the relative velocity of the wall with respect to the free-interface changes its sign. The 
transition in this part of the phase-diagram is continuous and belongs to the MN universality 
class. For q > q* the effective attractive force binds the interface to the wall and a continuous DP 
phase-transition takes place at a critical value v'^{q) < 1/2. 

In order to estimate the role of long range interactions in the critical dynamics of the SSW 
model, we have measured also in this case the effective activation rate X{£). Since the time interval 
1/vw between two consecutive wall moves can be regarded as the natural time scale of the SSW 
model, we measured the pinning rate at time tm = fn/vw (m = 1,2,...) for an initially depinned 
site located at the border of a depinned island of size H. at time tm-i • Averages have been taken over 
time, space and different ensemble realizations (typically 100) for large enough systems (L = 10^) 
close to the DP critical line, at g = 0.7 and q = 0.8. Our results are shown in Fig. 1111 where 
numerical data has been fitted with A(^) = A(l + a/£°"). There is evidence that A(£) decays to 
a constant faster than 1/^ with exponents a = 1.20(5) (for q = 0.7) and a = 1.30(5) {q = 0.8). 
According to our predictions based on the behavior of the generalized contact process, this implies 
a continuous DP transition, which indeed has been observed in Ref. ^. 

Unlike the RSOS model discussed in the previous section, here in the SSW model the DP 
behavior is observed at small length scales, and is thus accessible in numerical studies of finite 
systems. The reason is that the SSW model is designed in a way that the nonlinearity is maximal 
|32| . and hence the crossover length Ic is fixed and of order 1. On the other hand in the RSOS 
model the nonlinearity depends on the growth parameter with a diverging crossover length as 
p — > 1. 



V. CONCLUSIONS 



In this work we examined a possible connection between wetting phenomena and contact pro- 
cesses. Inspired by the puzzling richness of the phase-diagrams found in various models, we 
introduced a generalized contact process with the goal of capturing both the DP and first-order 
transitions observed in one-dimensional non-equilibrium wetting transitions occurring at non-zero 
interface velocity. The element of novelty distinguishing our model from standard contact processes 
consists in an algebraic dependence of the activation rate on the length i of the depinned island 
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containing the inactive site A(^) ~ A(l + a/£""). 

A mean-field analysis predicts that when cr > 1, the model exhibits a continuous phase- 
transition characterized by DP critical exponents. In other words, the algebraic decay of the 
interactions is not so long-range as to alter DP critical properties. Conversely, for < cr < 1 the 
phase-transition turns to a first-order one. Numerical simulations confirm the mean-field predic- 
tions, which are found to provide also an accurate estimate for the critical value CTc, separating the 
two different regimes: Cc = 1.0 it 0.1, from numerics. By directly measuring the effective activation 
rate A(£) in non-equilibrium wetting processes, one can use the estimate of the exponent cr as a 
practical tool for probing the nature of a wetting transition. 

To apply the insight gained from the DP process to wetting phenomena, we considered in this 
paper two wetting models, RSOS and SSW (see Section 4). Previous numerical studies of the 
latter model ^ indicate that for a sufficiently strong attractive force of the substrate (wah), the 
wetting transition is DP. Our criterion confirms these results, since the effective activation rate is 
indeed found to converge faster than l/£ to its asymptotic value. The phase-diagram of the RSOS 
model on the other hand is known to be more complicated, since it contains both first order and 
continuous wetting transitions 3]. Our analysis suggests that for p < 1 the entire phase-transition 
line located at the upper border of the coexistence region should asymptotically belong to the DP 
universality class. However the finite-size behavior depends on two length scales: (z) the crossover 
ic between EW and KPZ roughening behavior; (ii) the effective size Iq of depinned islands. As 
long as the system size is smaller then Lc = ic exp(^c/^o) and islands of size larger than £c are not 
significatively generated, numerical simulations can only provide evidence of a first-order transition. 
A qualitative phase-diagram for the RSOS model in the case qo < q^ip) is sketched in Fig. El 

It is interesting to compare our results with those recently obtained for the KPZ equation with 
an attracting hard core potential h = DV^h + v{yKf - V'{h) + C- Although the authors 
do not exclude the possibility that the entire transition line belongs to the DP universality class, 
their numerical analysis revealed both a first order and a DP wetting transition. In particular, 
DP critical behavior has been observed upon decreasing the attractive force and reducing by a 
factor ten both the diffusion (D) and the nonlinear (u) term with respect to the parameter values 
corresponding to a first-order phase-transition. While decreasing the attractive force is obviously 
expected to increase the "cut-off" scale io, the latter change decreases the length scale ic by the 
same factor ten. The crossover length from EW to KPZ roughening is in fact known to scale as 



33| . where T is the amplitude of the Gaussian white noise C- This indeed suggests that 



22 




0.5 



1.5 



FIG. 12: Qualitative "finite-size" phase-diagram of the RSOS model in the case qq < q^ip)- We expect the 
entire transition line for p < 1 to belong to the DP universality class. The full line (line A) on the left of the 
equilibrium point (EQ) marks a region in which the transient first-order behavior is the only numerically 
accessible one, with a crossover lenght Lc to DP behavior which diverges as, p 1. On the other hand, 
also systems of moderate size can exhibit a DP critical behavior over the dashed line (B). The dot dashed 
transition line (C) on the right of the equilibrium point is genuinely first-order (see text). Finally, the light 
dashed line where the free interface velocity changes sign marks the lower boundary of the phase-coexistence 
region (PC). 



the DP (first-order) behavior has been observed in a parameter range where the crossover scale 
Lc (as defined by Eq. ((T^ ) is numerically accessible (unaccessible), as suggested in the present 
study. This work helps to better understand the part of the wetting phase diagram characterized 
by a negative coefficient of the KPZ nonlinearity, which is precisely the part related with the 
complete synchronization transition (in the case of an MN phase transition, this relation can be 
made explicit by use of the Cole-Hopf transformation 3^). In particular our analysis indicates 
that the dynamical details of the system may induce a seemingly first order phase transition which 
lasts over exponentially long time and space scales, effectively suppressing DP critical properties. 



Finally, we wish to comment about the transition line for p > 1 (where the coefficient of the KPZ 



nonlinearity is positive). Forgo < ll 
the interface velocity changes sign 



)), the transition is known to be first-order, but it occurs when 
I, i.e. there is no region of phase-coexistence. Accordingly, 
depinned islands are no longer characterized by a triangular shape and interface fluctuations may 
easily give rise to the pinning of inactive sites far away from the active ones; in other words, the 
analogy with contact process is seemingly lost. It would be interesting to investigate whether the 
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inclusion of some sort of "spontaneous" nucleation of active sites can eventually account for the 
scenario observed for p > 1. A further open problem is the crossover from the DP to the MN 
universality class, which takes place in both RSOS and SSW models. The study of this crossover 
would require the inclusion of nucleation of active sites in the interior of inactive domains. An 

appropriate generalization of the DP model including such processes could yield useful insight onto 
the crossover phenomena taking place in non-equilibrium wetting. 
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